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ABSTRACT 

We revisit the calc ulation of the Ohmic dissipation in a hot Jupiter presented 



m 



Laine et al.l ( 120081 ) by considering more realistic interior structures, stellar 
obliquity, and the resulting orbital evolution. In this simplified approach, the 
young hot Jupiter of one Jupiter mass is modelled as a diamagnetic sphere with 
a finite resistivity, orbiting across tilted stellar magnetic dipole fields in vacuum. 
Since the induced Ohmic dissipation occurs mostly near the planet's surface, we 
find that the dissipation is unable to significantly expand the young hot Jupiter. 
Nevertheless, the planet inside a small co-rotation orbital radius can undergo 
orbital decay by the dissipation torque and finally overfill its Roche lobe during 
the T Tauri star phase. The stellar obliquity can evolve significantly if the mag- 
netic dipole is parallel/anti-parallel to the stellar spin. Our results are validated 
by the general torque-dissipation r elation in the prese nce of the stellar obliquity. 



We also run the fiducial model in lLaine et al.l (120081 ) and find that the planet's 



radius is sustained at a nearly constant value by the Ohmic heating, rather than 
being thermally expanded to the Roche radius as suggested by the authors. 



Introduction 



In the study of the orbital distribution of known Jupiter-mass exoplanets, the radial- 
velocity method has revealed a pile-up of hot Jupiters with orbital periods of ~ 3 days 
(e.g., see The Extrasolar Planets Encyclopedia website at |http://exoplanet.eu[ ). A number 
of models have been propo sed to explain t he pile-up, such as an inner disk cavity stop - 



ping planet migration (e.g. iLin et al.lll996t iRice et al.ll2008t iBemtez-Llambay et al.ll201ll ). 
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and tidal circularization of a gas giant planet in an extremely eccentric o rbit arising from 



planet-planet interactions after the proto-plane t ary disk dispers es (e.g. IChatterjee et al. 



20081 : iNagasawa et al.ll2008l : IWu fc Lithwick 



20 lit iNaoz et al.ll201lh . In addition, tidal heat- 



ing in a young hot Jupiter in a moderately eccentric orbit may inflate the planet over its 



Gu et al. 


2003, 


2004; 


Chang et al. 


2010) 



tation of the planet's ec centricity in this case is subject to the uncert ain density profile of 



the inner edge of a disk ( iRice et al.ll2008t iBemtez-Llambay et al.l 120111 ). 



Aside from t he tidal dissipation that relies on the presence of orbital eccentricity, 
Laine et al.l ( 120081 ) invoked the Ohmic dissipation to inflate a young hot Jupiter in a circular 
orbit by adopting the model proposed by Campbell (1983,1997) for the magnetic interactions 
in the AM Herculis systems. In this simplified model, the planet is assumed to behave as an 
imperfect conductor without its own ionosphere and magnetosphere; namely, a diamagnetic 
sphere with a finite resistivity. In addition, it is assumed that the stellar spin is aligned with 
the planet's orbit. Since the vacuum space is assumed between the star and the planet, the 
stellar magnetic dipole must be misaligned to induce electric currents and magnetic fields as 
the planet circles its T Tauri star. The magnetic torque arises from the Ohmic dissipation 
in the planet at the expense of the spin-orbit energy (see the §2.3). 

Normally a planet even without its own fields possesses an ionosphere due to the photo- 
ionizati on of the upper at mosphere. An induced magnetosphere can form above the iono- 
sphere (jZhang et al.l 120091 ) . As a planet orbits a star with a tilted magnetic dipole, the 



ionosphere may shield the time-varying stellar fields so sufficiently that little electromag- 
netic field can be induced in the planet's interior by the external stellar fields. Nevertheless, 
one may argue that a young hot Jupiter is already tidally locked by its parent star such that 



circulation in the atmosphere (e.g. 


Showman et al. 


20081. 


2009 


• 


Thrast arson & Cho 


2010; 


Rauscher & Menou 


2010; 


Dobbs-Dixon et al. 


201(1 1 


^erna et al. 


2010), which may maintain 



an ionosphere on the permanent night side. Based on the radio-sounding results from the 
Venus Express spacecraft, the ionosphere on the night side of Venus is weaker and possibly 
more sporadic than that on the day side (IPatzold et al.l 120071 ). Hereafter, we boldly apply 
the model for AM Her binaries to the entire planet and also consider a vacuum space outside 
the planet and the star to simplify the calculation. The consequ ence of this simplification 
is that other electromag netic effects such as unipolar induction (iGoldreich fc Lvnden-Bell 
19691 ; iLaine fc Linl 120121). Alfven- wave w ings (|Neubauerlll980l ; iKopp et al.l 120111). dynamical 
friction (Papaloizou 2007 ). stellar winds ( Vidotto et al. 2010 ). planetar y winds dAdams 2011 



Trammell et al 



2011 ). stellar fields diffusing into the planet interior ( Campbell 20051 ). and 



magnetic reconnections are all ignored (also see Lanza 2011 for a recent review). In addition, 



-3 - 



any electromagnetic effect s associated with atmospheric circulations are not being consid- 



ered for further simpli city ( IPerna et al.ll2010l ; iBatygin fc Stevensonll2010l ; iBatygin et al.ll2011 
Wu fc Lithwickl |2012| ) . In short, we restrict ourselve s only to the diam agnetic part of the 
star-planet magnetic interaction^, as was modelled bvlLaine et al.l (120081 ) . It should be noted 



that the Ohmic heating proposed by lLaine et all J2008h is short-lived, since the stellar mag- 
netic fields decay sig nificantly during the T Tauri phase. This is in contrast with the Ohmic 
dissipation model by IBatygin et al.l ( 1201 if ), which is long-lived. Consequently, this study is 
concerned exclusively with the early evolution of hot- Jupiter systems. 



The observations using the Rossiter-McLaughlin effect (jOhta et al.l 120051 ) suggest that 
dwar f stars hosting transiting planets may have possessed a wide range of stellar obliquities 



^e.g. iWinn et al.l |2010| . 1201 ll ). These observational findings seem against the conventional 
paradigm in which a planet should orbit in the same direction as the stellar spin as the star 
and planets form together in a proto-planetary disk. A number of N-body numerical simula- 
tions demonstrated that after the proto-planetary disk disperses, planet-planet interactions 
accompanied by tidal circularization, as mentioned in the first paragraph of the Introduc- 
tion, can generate obliquities. It was also proposed that before the proto-planetary disk 
dissipates, the warp torque resulting from the magnetic interactions between the proto-star 
and the inner part of the disk would move the stellar spin away fr om the disk angular mo- 
ment um despite the presence of gas accretion onto the proto-star ( ILail 12012c iFoucart fc Lai 
201ll ). Motivated by the latter works, it is timely to consider a more complex case in which 
stellar obliquity A is not zero; i.e., the orbital axis is not aligned with the stellar spin. 

It should be noted that the tidal dissipati on in the star drives the system to the spin - 
orbit alignment as well as synchronization (e.g. lHutlll98ll ; iMatsumura et al.ll2010l ; lLaill2012l ). 
To make the problem tractable, we do not take account of the influence on A driven by 
the proto-planetary disk or by tidal interactions with the proto-star, but simply take A as 
a free parameter in this work. In addition, we assume that the planet spin is tightly being 
synchronized with its orbital motion during the evolution, therefore generating negligible 



dissipation in the planet (e.g. iGu et al.l 120031 ). This simplification allows us to ignore the 



effect due to the planet spin in the calculation. 

Owing to the Ohmic dissipation and the resulting magnetic torques, the stellar spin, 
planet's orbit, and the interior structure of the planet evolve simultaneously. To calculate 
the coupled evolution more precisely, we adopt an interior-structure model (for details see 



1 The same concept has been applied to star-disk magnetic interactions 
response of the disk is modelled by a diamagnetic disk as well as a magnetically threaded disk. Unlike our 
model planet possessing a finite resistivity, the disk is assumed to be a perfect conductor in the diamagnetic 
part of their model. 



Laj|[l999) in which the magnetic 
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Chang et al.ll201Q[ ) to compute the planet resistivity and the thermal response of the planet 



due to the Ohmic heating. 

The structure of the paper is organized as follows. In §2, we describe the equations for 
the coupled secular evolutions of stellar spin, planet's orbit, and planet interior structure 
due to the diamagnetic interaction between a young hot Jupiter and its parent T Tauri 
star. In order to understand the dependence of Ohmic dissipation on various orientations of 
the stellar spin and magnetic dipole moment, we first conduct a parameter study in §3 to 
investigate this with no secular evolutions. The parameter study involving secular evolutions 
is then presented in §4. Finally, we summarize and discuss the results in §5. 



2. Governing equations for the coupled evolutions of spin, orbit, and planet's 

interior structure 



Following the same mathematical procedures in lLaine et al.l ( 120081 ). we solve the resistive 
induction equation in the co- moving frame of the planet, with the stellar dipole fields and 
the induced fields expressed in terms of the poloidal scalars 0*(r, t) and 4> p (r, t), respectively: 
namely, the magnetic field B, which has a poloidal nature in our problem, is related to the 
poloidal scalars by B = V x (V x (0f)), where r is the unit vector of r. The SI unit system 
is adopted to present the equations for electromagnetic calculations. The induced poloidal 
scalar p can be solved by the separation of variables in the spherical coordinates (r, 9, if) of 
such a frame after 0* and the resistivity profile rj(r) are given. Let R p be the planet radius. 
For notation convenience, we denote p (r > R p ) = 4> P)OU t and p (r < R p ) = 0j n . Hence the 
total poloidal scalar outside of the planet is <p out = 0* + 4> P)0 ut- In the case of A = 0, 0* and 
therefore P vary at the rate equal to w t - n as viewed by the planet, where is the stellar 
spin angular frequency and n is the orbital angular frequency of the planet. In the Appendix 
A, we illustrate the coordinate systems for the problem (see Figured]) and derive the detailed 
equations to solve for the potential scalar induced by a tilted magnetic dipole in the presence 
of stellar obliquity. We show that in order to describe the time-varying potential, there will 
be 3 more frequencies involved other than w_ = us* — n; they are co + = + n, and n. 
Once the potential scalar P is solved, the induced magnetic field = V x (V x p f), the 
electric field E, th e electric current j, and hence the Ohmic dissipation can be all calculated 



(lLaine et al.ll2008l ). In the following subsections, we describe how to calculate the resistivity 
of the planet and the corresponding spin and orbital evolutions due to the Ohmic dissipation 
in our model. 
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2.1. Calculation of resistivity 



For a hot Jupiter around a T Tauri star of one solar luminosity, the equilibrium tem- 
perature is ~ 1000 — 2000 K at the photosphere. The gas in the region just below the 
photosphere is therefore weakly ionized due to thermal ionization of alkaline elements. As 
the temperature and density continue to rise in deep layers, the thermal ionization of the 
most abundant constituents H and He starts to become non-negligible. In the even deeper 
interior, the density is hig h enough so that the fluid is partially degenerate and fully ionized 
due to pressure ionization (ISaumon et al.lll995l ). It has been shown that the electric currents 
and magnetic fields induced near the planet's surface are only present in the outer part of 
the planet where the ionization fraction is low and hence the resistivity is high. That is, 
magnetic fields decrease significantly oyer a skin depth S from the surface to the interior 
(ILaine et al.ll2008l ; iBatygin fc Stevensonll2010l ). In other words, the induced electric currents 
and magnetic fields are considerably shielded out by the outer part of the planet such that 
the precise values of resistivity in the interior do not matter. Thus in this work, we restrict 
ourselves to the resistivity 77 due to elect ron-neutral collisions in a weakly ionized plasma 
(braine et al1ll983l : iBlaes & Balbulll994h : 



Ve . n = 230 fcf) T 1 / 2 cm 2 /s, 



(1) 



and apply the above equation to the entire planet without making a significant error. In the 
above equation, n n is the neutral number density, n e is the electron number density, and T 
is the temperature. 

To estimate the ionization fraction in eq.([T]), we first consider the ther mal ionization of 
alkal i ne elements. Ther mal ionization is governed by the Saha equation (cf. IBlaes &: Balbus 
19941 : IPerna et al.lhoioh 

m 3/4 



(2) 



where n e = Yl n ji n j = fj n i k * s the Boltzman n's constant, h is the Planck 's constant 
divided by 2n, and <C rij is assumed. We follow IBatygin fc Stevenson! (120 lOl ) to find the 
abun dances and ioniza tion potential of each a lkaline species (labelled by j) inferred 
from lLoddersI J1999H 2 ! and ICox fc Pilachowskil feoooh . 



2 fj are estimated at the temperature > 2000 K. Below this temperature, the abundances of some species 
such as Fe and Ca decline dramatically due to their molecular formations with other atoms. This process 
does not affect our results significantly because K and Na are the primary sources of thermal electrons at 
the low temperatures. 



-6 - 



In even deeper layers, the thermal ionization of the most abundant constituents H and 
He starts to dominate the electro n contribution . We c ompute the H & He ionization based 
on the equation of state tables in lSaumon et al.l ( 119951 ). Given th e pressure P, tempe rature 
T, and the helium mass fraction Y, the mass density is given by (jSaumon et al.lll995f ) 



1 



1 - Y 



Y 



p(P,T) p H (P,T) p He (P,TY 



(3) 



where p H {P,T) and p He {P,T) are obtained from interpolation of the data in the EoS tables 
for pure H and He, respectively. Hence, the electron number density n e and the total number 
density n are given by (see eqs.(36) & (37) in Saumon et al. 1995) 



n 



2p H /rriH _ x h + 



1 + 3X H2 + X 



3pHe/m He xIIe 

1 + 2Xff e + X}j e + e 



f^H 2 ,H,H+,e + nHe,He+,He 2 +,e 



2p H /m H 



3Xff 2 + X 



+ 



3p He /m He 



H 



2X He + X 



He+ 



(4) 
(5) 



where pn ~ p(l — Y), pue. = P — Pn = Yp, X^ and X^ e are given by eqs.(34) and 
(35) respectively in ISaumon et al.l ( 119951 ). Y = 0.283 is adopted in our interior-structure 
simulations. 



2.2. Planet radius and spin-orbital evolution due to Ohmic heating 



Including the Ohmic dissipation but neglecting the small planetary spin energy (IBodenheimer et al 



200 i f), we have th e evolution of the global energy for the entire planet governed by (cf. 
Chang et allboioh 

U + W = Q ohmic - L, (6) 

where U is the internal energy, W is the gravitational potential energy, L is the intrinsic 
luminosit y from the photo sphere of the planet, and Qohmic is the Ohmic dissipation rate 
given by (jLaine et al.l 120081 ) 



Q 



ohmic 



p oV Re(i) 2 dV 



>r<R v 



[Re(V x B)] 2 dV 



(7) 



where "Re" means taking the real part and () denotes the time averaging over the time scale 
longer than the forcing periods; namely, it is the secular evolutions of the spin and orbit that 
are relevant to the long-term thermal evolution of the interior structure. Qohmic should be 
equal to the average flow of the electromagn etic power (i.e . the Poynting vector) into the 
planet through the planet's photosphere (e.g. |Jacksonlll990l ). 
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Owing to the diffusive nature of the problem, an order-of-magnitude estimate for the 
Ohmic dissipati on can be made based on the dimen sional analysis of Equation[7|with V ~ 1/8 



and dV ~ E*6 f lCampbeimi997t lLaine et all 120081 ) 



B 2 

Qohmic ~ — (AttR 2 5)uj, (8) 
2/x 

where the skin depth for our magnetic induction problem is 5 = (277/w) 1 / 2 , the stellar 
magnetic field near the planet is B ~ (/io/47r)m/-D 3 , u is the forcing frequency, m is the 
magnitude of the stellar dipole moment, and D the orbital separation (see Figured]). 

The spin-orbit evolution is dictated by the dissipation torque. The torque acting on the 
stellar spin due to the electromagnetic interaction in the inertial frame takes the form (cf. 



Campbell! 119971 ) 



Tinert l^obliquity, inert X Bp] ane ^(r I" s ^ ar ) m o h,j;q U ity , inert X V (*9r0p)r=r s tar • (^) 

In the above equation, m ob i iquity iner t, given by P x 'Hi (see the Appendix A), is the stellar 
dipole moment as seen in the inertial frame with n in the ^-direction, and B p i anet is the 
planet-induced magnetic field at the location of the star r star — Di — Dr + (ir/2)0 (see 
Figure 1). Note that although B p i anet is calculated in the planet's rest frame, its value 
in the inertial frame is the same as in the non-relativistic regi me when the terms of order 



n 2 D 2 /c 2 <C 1 are neglected, where c is the speed of light (e.g. iThyagaraja fc McClements 



2009|). 



We shall see that the stellar spin precesses secularly with time in the inertial frame 
and we are interested in the spin-orbit evolution on the secular timescale, much larger than 
the spin and orbital periods. This amounts to taking the time-average for each physical 
quantity to average out their short-term variations. Besides, it is convenient to work out 
the secular evolution problem in the precession coordinates with n always pointing to the 
z-direction. Hence after each time step, we switch to the inertial frame such that the stellar 
spin is always on the y'-z' plane and the orbital angular momentum is always along the 
z'-axis at the beginning of the next time step. We denote this "instantaneous" inertial frame 
as 0'x"y"z', which coincides with O'x'y'z' at the beginning of e ach time step. Thus the 



time-averaged torque (T; nert ) in this inertial frame is given by (cf. ICampbelllll997l ) 



\T X ") = m{rh z 'B p i anettr — rhyi'B p i ane t t 0) , (10) 

(Tyif) = m{rh x ii Bpi ane t + m Z ' Bpianet^), (11) 

\T Z ') = —m(m x iiBpi anet , r + rhy>'B p i anettlfi ), (12) 
where m x ", m y », and m z i are the three Cartesian components of the unit vector of m ob i iqu i tyjinert . 
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Once the time-averaged torque is obtained in the inertial frame, we are ready to calculate 
the secular evolutions of spin and orbit. First of all, (T x ii) leads to the precession of stellar 
spin around the orbital axis; namely, 

d(<p") 



/(A)J*o;*- 



dt 



(T x »), 



(13) 



where (ip") refers to the time-averaged precession angle and /(A) is a function of A given 
in the Appendix B. In the "instantaneous" inertial frame, (Ti ner t) along the stellar axis 
determines w*. — (T z i) governs h and thus D. Moreover, A is caused by the components 
of (Ty/i) and (T z i) normal to the stellar spin, and additionally by the back reaction — (T y ») 
acting to the orbital angular momentum. In other words, using the "instantaneous" stellar 
spin uj* = w*a>* = co>*(sin Aj" + cosAk ), we arrive at a set of evolutionary equations: 

d(I*uj* 



dt 



d\ 
~dl 

M, 



(Tinert) ■ = (T y ») sin A + (T z >) COS A, 
(T y rr) cos A — (T z i) sin A (T y ") 



d(D 2 n) 
' dt 



M p D 2 n 



(14) 
(15) 

(16) 



where M p is the planet's mass and the complicated expression for dX/dt is explained in 
the Appendix B. The above 3 equations can be combined to express A in terms of L spin = 
d(Lu*) / dt and L or b = M p d(D 2 n) /dt as follows 



d\ 
~dt 



dL. 



orb 



dt 



Leu* sin A 



+ 



M p D 2 n tan A 



+ 



dL 



spin 



dt 



1 



Lu* tan A 



+ 



1 



M p D 2 n sin A 



(17) 



which we shall find quite useful to interpret the evolutionary results. 



Note that the moment of the inertia L of the T Tauri star also evolves. In reality, m 
evolves as well (j Johns- Krulll 120071 : lYang fc Johns-Krullll201ll ). but in this work we prescribe 
a constant value for m for simplicity. 



When A = 0, the terms in the magnetic potential scalar 
are left, leading to (T x n) = (T y ») = 0. Therefore Equations (CC 



associated only with w_ 
-fTTBT) reduce to the ones in 



Campbell! fll983f ): 



M, 



d(Luj*) 

dT~ 

d(D 2 n) 

' dt 



~(T zl ). 



(19) 



In the absence of stellar obliquity, the torque and Ohmic heating can be simply related 



to each other by virtue of the equation Q 



ohmic 



\u_(T z /)\ (jCampbelll Il983t lLaine et al. 
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20081 ). If A 7^ 0, we show in the §2.3 that the vector product (Ti ner t) • (<^* — n ) gives rise to 
the Ohmic dissipation that drives the spin-orbit system toward a lower energy state. More 
specifically, 

Qohmic = ~ {u*(T y ») sin A + u*(T z >) cos A - n(T z >)) , (20) 

= - U *^r +n ^r- (21) 

Since Qohmic, w*, and n are all positive quantities in this work, the above equation indicates 
that dL spin /dt > dL^/dt for u;* < n and vice versa, which is a familiar result for A = 
but even applies generally to the cases for A ^ 0. Note that even when n = a;*, Qohmic ^ 
due to the spin-orbit misalignment. In addition, (T x >t) causing the precession of stellar spin 
around the orbital axis does not do any mechanical work and thus is not related to the 
Ohmic dissipation in the planet. 

In the special case where a = 0, the stellar dipole moment in the inertial frame is 
rti biiquit y ,inert = (0, sin A, cos A) . Substituting this into eqs. flTT]) & (Tl2|) . we have the unique 
relation (T z r) / (T y ") = — tan A regardless of the value of uj*. This together with eq. (!20|) gives 

Qohmic = n(T z i), (22) 

which is independent of as it should be when the spin and stellar dipole are aligned. Once 
again, Qohmic > and n > by our sign convention. It follows from the above equation that 
(T z >) > 0, therefore always leading to an orbit decay. 



2.3. General relation between energy dissipation and torques 

Since the torques arise from energy dissipation, we wish to derive the relation between 
the torques and dissipation in the presence of obliquity. This relation provides a powerful 
check on whether the Ohmic dissipation and the resulting torques calculated in §2.2 are 
correct. 

The stellar spin angular momentum and planet's orbital orbital angular moment are 
given by 

L S pin — (23) 

and 

Lorb = {m p \J Gm*a)ri, (24) 
respectively. Since the total angular momentum is conserved, T" = dh sp i n /dt = —dL or b/dt. 



-lo- 



pe) 



However the total energy of the system is not conserved as a result of dissipation. The 
stellar spin energy changes at a rate according to 

dE spin _ d(l/2)I*u>l _ ^ dI*co* _ ^ dL spin 
dt dt dt dt 

The change rate of the orbital energy is 

dE orb GM^nip da 
dt ~ 2a 2 dt' 

which can be linked to the change of the orbital angular momentum as follows 

dE or b _ <iL or b ^ 
dt dt 

In deriving the above equation, we have taken the time derivative of Equation (T24|) and used 
the identity dh/dt ■ n = 0. Thus 

• dyE S pi n -\- Eq^J DLspin . ^Lorb m// / \ / nn \ 

-Q ohmic = Jt = ^. — + — .n = T >,-n). (28) 

Although the heating rate is expressed in terms of the Ohmic dissipation, the above relation 
for dissipative torques can apply generally to other dissipative processes such as tidal dis- 
sipation. It is apparent that we do not specify the Ohmic dissipation to deduce the above 
relation. 



2.4. Summary of procedures 

Given E(t) and M p , the evolutions of stellar spin and orbit (a;*, A, and D) are coupled 
with the evolution of interior structure (R p , rj(r) etc.) via the Ohmic dissipation in a hot 
Jupiter, which is modelled as a diamagnetic sphere in our calculation. 

The procedure of the evolutionary calculations is summarized as follows. We start 
with initial a;*, A, D, and R p to obtain Qohmic and T inert . The next step consists of 
three calculations: the first is the calculation of the new interior structure of the young hot 
Jupiter due to Qohmic, the second is the computation of the new E from a stellar code, 
and the last is the calculation of the new w*, A, D from the integration of the ODEs from 
Equation (TH]) to Equation (ITS]) based on T inert . Consequently, i*, w*, A, D, and R p at 
the next time step will be obtained. Meanwhile, the computed T inert and Qohmic can be 
checked using Equation (120]) to validate the calculation. The same procedure is then carried 
out over and over again to evolve the system until either the planet reaches its Roche radius 
or the calculation approaches the end of simulation at 10 7 years. We employ the same codes 
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described and used by iBodenheimer et al.l ( 120011 ) and IChang et al.l (120101 ) for the planetary 
and stellar interior structures, respectively. 

To simulate t he stellar rotation being locked by a process such as disk locking (e.g., see 
Chang et al.ll2010l and reference therein), we also run cases (actually most of the cases) in 
which the stellar spin is held at its initial value throughout the simulation, even though 
A is still allowed to evolve by the magnetic interaction. It is conceivable that any external 
torques af fecti ng uj* should change A as well, such as the star-disk magnetic interaction by 
Lail (120121 ) and lFoucart fc Lail (120111 ). In this study, the star-disk interaction is not modelled 
with the star-planet magnetic interaction. Instead, we focus only on the evolutions due to 
the star-planet magnetic interaction, with the condition for to be "locked" for the sake 
of simplicity of the toy model. 



3. Comparative studies without secular evolutions 

In this section, we present a couple of test runs in our model without considering the 
evolution of spin, orbit, and interior structures of the proto-star and planet. The purpose of 
the test runs is to investigate how the Ohmic dissipation varies with A and a. This provides 
parameter and thus comparative studies to understand the basic behavior of the results 
before we proceed to the more complicated calculations involving secular evolutions. 

Such comparative studies for the Ohmic dissipation rate Q ohmic vs. uo*/n are shown in 
Figure El Q ohmic is calculated by virtue of Equa tion (J7J). We adopt m = 4 x 10 34 A m 2 , 



the same fiducial value used in lLaine et al.l (120081 ) . The interior structure for a coreless hot 



Jupiter with M p = Mj and R p = 1.84Rj is used for the test runs. We consider as a free 
parameter, whereas n is held constant corresponding to the orbital radius of 0.02 AU. Thus, 
the stellar irradiation, which affects the interior structure, is also constant. This reduces the 
number of variables and helps to more easily examine how Qohmic varies with u*/n in the 
test runs. Everything else being the same, there is no difference in the Ohmic dissipation 
for a and for (180° — a) due to the axi-symmetry of dipole fields. Consequently, we only 
present the cases for a < 90° in Figure El 

The upper left panel of Figure El shows that in the absence of the stellar obliquity 
(A = 0), the Ohmic heating rates increase from zero for a = to the maximum values for 
a = 90°. In addition, the heating rate vanishes when w_ = (i.e. u*/n = 1) and increases 
with due to stronger electromagnetic interactions induced by f aster forcing. The trend 



and the character of these results agree with those in lLaine et al.l (120081 ). The similar line 



of argument applies to the cases for A = 180° in which the stellar spin is completely flipped 
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over and therefore the forcing frequency is uj + rather than As illustrated in the lower 
right panel of Figure [2], the Ohmic heating rate increases with a. Besides, the heating rate 
increases with and thus u + . 

We find that the dissipation torque \T Z >\ for A = 0° decreases with the forcing frequency 
\ui-\ except when the forcing frequency is very close to zero; i.e. the torque peaks at w_ ~ 
1.3 x 10 -4 s -1 for A = 0°. The maximum value of the torque arises because the torque is 
significantly weak for extremely slow forcing, and becomes small again for fast forc ing due to 



the dissipa tion localized within one small skin depth below the planet's surface (jCampbell 



19831 . Il997h . 



When A 7^ 0, the relation between and Qohmic becomes perplexing and requires more 
explanations. As shown in Figure [2] for A = 45°, 90°, and 100°, the positive correlation 
between and Qohmic exists when u*/n is large enough for the forcing frequency to play 
the main role. This outcome can be realized by contemplating the problem in the two extreme 
regimes: a > and a < 90°; the orbital motion alone contributes most of the heating in the 
former regime, whereas in the latter regime the Ohmic heating is generated primarily from 
the relative spin-orbit motion (i.e. w_ or u + depending on A). More specifically, Figured 
shows that the heating rate is constant independent of for the cases of a = 0, in agreement 
with Equation fl22|) ; namely, the Ohmic dissipation induced entirely by the orbital motion 
with the forcing frequency n. As a starts to deviate from zero, we find that the Ohmic 
dissipations induced by other forcing frequencies begin to increase but the Ohmic heating 
arising solely from the orbital motion starts to decrease. This can been seen in Figures [3] & 
H]for a = 10° where the total heating profile (shown in cyan line) near the planet surface 
almost overlaps with the one corresponding to the forcing frequency n (dotted blue line), 
but the heat contributions from other frequencies other than n are not totally negligible. 
When a = 90°, the heat contribution from the forcing frequency u + or w_ as a result of the 
relative spin-orbit motion totally dominates over that from the forcing frequency n; namely, 
in the outer part of the planet, the total heating profile (cyan line) almost coincides with the 
one for U- (dashed green line) in Figure |3] and for u + (magenta line) in Figure HI It can be 
confirmed by Equations f lAip and (|A3[) that when a = 90°, the forcing with the frequency 
n disappears in the expression of the stellar magnetic dipole moment m ob i iquity and thus in 
the corresponding poloidal scalar 0*, leading to null contribution of the dissipation from the 
forcing frequency n. Therefore, the tiny dissipations for a = 90° shown in Figures |3] & H] 
(i.e. dotted blue line) stem totally from numerical errors, which are too small to affect the 
results. Note that the dissipation occurs primarily in the outer part of the planet because the 
induced magnetic fields are mostly confined within one skin depth below the photosphere. 



Given the above explanations, we are able to further elaborate the general dependence 
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of Qohmic on A and a shown in Figure [2j Let's first examine the cases for A = 45°, 90° and 
100°. In these cases, the heating rate in the high frequency range u*/n > 1.5-2 increases 
with a. Roughly speaking, it is because the primary forcing switches from the slow rate n 
to the fast rate u + (for A > 90°) or uj_ (for A < 90°) as a increases (see the oj^/n = 2 case in 
Figure [3]). The trend apparently reverses in the low frequency range w*/n < 1.5-2; namely, 
the heating rate decreases with the increasing a (see the uj^/n = 0.5 case in Figure [3]). 

In contrast, for the retrograde orbits with large stellar obliquities as represented by the 
cases for A = 135° and 180° in Figure [2J the Ohmic dissipation always increases with a. The 
forcing with the frequency u + is always fast enough to induce more heat for larger a than 
the heat generated mostly by the slower forcing with the frequency n for smaller a. This 
consequence can be implied by comparing the heating profiles for u^/n = 2 with those for 
u*/n = 0.5 in Figure HI 

To further validate our numerical calculations, we also compute Qohmic based on the 
general torque-dissipation relation given by Equation ( 1201) and show the results in Figure |5j 
In general, Figure[2]and Figure[5]are consistent with each other. The discrepancy between the 
two types of calculations of Qohmic is < 7%. In addition, the results of Qohmic can be crudely 
verified by Equation ([8]). Using rj = 8.3 x 10 9 m 2 /s, which is approximately the maximum 
value of the 7](r) profile in the calculations, we obtain the skin depth 5 ~ 1.5 x 10 9 cm. The 
substitution of this skip deptbH into Equation flH]) gives the dissipation rate Qohmic ~ 2 x 10 31 
erg/ s, which is on the similar order of the magnitude to those shown in Figure |2J 

All the dissipation rates in the test runs have been calculated based on Equations ([I]) 
and fl2]) under the assumption of the low ionization fraction of each alkali species as well as 
the low total ionization fraction within the skin depth. To verify whether this assumption is 
reasonable in terms of the heat generation , we a pply the full version of the Saha equation, 



e.g. Equation(l) in iBatygin fc Stevenson! ( 120101 ) . to the test runs. We find that the total 



ionization fraction is sufficiently low in the outer part of the planet such that Equation ([I]) 
still applies. We then compute the new heating rate profiles and compare them to those 
based on Equation (j2]) . Figure [6] illustrates the comparisons for the two cases shown in the 
upper left and lower right panels of Figure El as the representative examples. It is evident 
from the figure that the heating profiles derived from Equation ([2]) and from the full version 
of the Saha equation are almost the same in the outer part of the planet where most of the 
dissipation occurs. It then follows that the total heating rates derived from the full version of 



3 In these calculations, the radiative-convective interface lies at about 1.295 xlO 10 cm; i.e. 1.9 x 10 8 
cm below the photosphere. Hence, the main heating region, characterized by the 5, extends down to the 
convection zone. 
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the Saha equation are only 3-4% higher that those derived from Equation fl2]), thus validating 
the approximate results using Equation (j2J). 



4. Evolutionary results 

We now present the evolutionary results with the input parameters and different initial 
conditions listed in Table [TJ The initial is given by one half of the initial n. This initial 
condition is based on the assumptions that the inner edge of the disk is loc ated at the location 



of the co-rotation radius of the proto-star due to disk locking (e.g., see iChang et all 120101 . 
and reference therein) and that the initial location of the planet lies in the orbit with the 
2:1 mean motion resonance with the inner edge of the disk according to planet migration 
theories (see Lin et al. 1996; Rice et al. 2008; cf. Bemtez-Llambay et al. 2011). The 
simulation is run from ti = 0.7 to t en d = 10 Myrs, corresponding to the T Tauri star phase. 
The starting time 0.7 Myrs is compara ble to the times cale of the type II migration time of 



a giant planet in a protoplanetary disk ( )Lin et al.lll996l ). Except for Case 1 which allows u 



to evolve according to Equation ffT4|) for comparison, we do not evolve in other cases as 
it is assumed to be locked by some process su ch as disk locking. We also run Case 20 with 



the parameters similar to the fiducial model in lLaine et all ( 20081 1: M p = 0.63M/, Di = 0.04 
AU, lu_ = 10~ 5 s -1 , A = 0°, a = 90°, M* = M , and L* = 1.5L . It should be stressed 
that even in an aligned system, the size of the magnetospheric inner cavity is proporti onal 



to m 4 / 7 M 2 / 7 , where M is the disk gas accretion rate onto the T Tauri star (e.g. see lLai 



19991 . and references therein). In other words, the initial n is in fact related to m and M. 
Moreover, evolves as the magnetospheric cavity evolves even in the disk-locking model. 
Since we assume a constant m and do not intend to model the cavity size in the presence 
of stellar obliquity and the misaligned magnetic dipole, we simply parameterize the initial 
value of n independent of m in this work. 

Cases 1-9 represent the evolutions of a young hot Jupiter of 1 Mj initially at the very 
close distance Di w 0.02 AU, resembling a planet lying inside a small magnetospheric cavity. 
Owing to the small starting orbital distance, the strong Ohmic dissipations are generated 
on the order of 10 30 ~ 31 erg/s throughout the evolutions, resulting in fast orbital decays. 
However, because the intense heating occurs mainly near the planet surface, the dissipation 
is unable to significantly inflate the planet against self-gravity. Figure [7] shows that the rise 
in R p is < 3% in these cases when the planet quickly shrinks its orbit and fills its Roche 
lobe in just a few 10 5 to about 1 million years after The small increase in R p arises from 
the thermal expansion of the outer part of the planet. Despite the intense heating near 
the planet surface, temperature inversion is not observed because the strong dissipation is 
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limited to the radiative layer in the evolutionary cases and thus is easily losto In these cases, 
the planet first undergoes relatively fast expansion as the dissipation is suddenly deposited 
in the beginning, and then reaches an intermediate quasi-equilibrium state (i.e. Qohmic ~ L 
in Equation [H]) that lasts for some period of time depending on how fast the orbit decays. 
The planet expands again as the orbit continues to shrink and thus the Ohmic dissipation 
is further enhanced. 

Among these cases, Cases 1,2, and 3 present the evolutions in the absence of the stellar 
obliquity (A = 0). In Case 1, evolves, caused by the dissipation torques and /* according 
to Equation (fT4"|) . without the spin-locking assumption. In Cases 2 and 3, cu* is constant. 
Figure [7] shows that R p and D in Case 1 evolve faster than those in Case 2, starting from the 
same initial conditions. It is because the forcing frequency \oj_ \ is lower in Case 1, leading to 
faster orbital decay as explained in §3. A larger skin depth results from the slower forcing, 
generating deeper heating and thus faster expansion. Figure [7] also shows that R p and \D\ 
in Case 3 are always larger than that in Case 2 as expected from the test runs in §3; the 
larger a in Case 3 produces the stronger heating and faster orbital decay. 

On the other hand, Cases 4, 5, and 6 present the studies in which A ^ but the stellar 
spin and dipole are parallel (i.e. a = 0). Although the stellar spins in Cases 4 and 6 point to 
opposite directions, Figure [7J shows that their evolutions of R p and D are similar due to the 
similar time variation of the stellar dipole field that is axi-symmetric about the spin axis. 
Furthermore, the larger R p and \D\ in Case 5 than those in Cases 4 and 6 is a result of the 
larger A and thus stronger heating, in accordance with the results of the test runs shown in 
the upper middle and upper right panels of Figure [2] for a = 0° and uj*/n = 0.5. The almost 
symmetric evolutions between Cases 4 and 6 are broken when a ^ 0°, as illustrated by the 
different evolution curves for their counterpart cases 7 and 8. 

We also run Case 9 to compare with Cases 4 and 7 to examine the evolutions starting 
from the same Aj = 45° but different a. As has been demonstrated in §3, Case 9 lies in 
the special regime where a = 90° and hence no dissipation is contributed from the forcing 
frequency n, in contrast to the other extreme regime shown in Case 4 where the dissipation 
is totally from the forcing frequency n. The total dissipation rate and therefore R p as well as 
\D\ increases with a in these cases during the evolutions, which is consistent with the result 
of the test run for A = 45° and uj*/n = 0.5 displayed in Figure [2j 



4 It is different from the non-evolutionary test runs shown in Figures [3] and [4] where the heating profiles 
extend down to the convection zone. When allowing for evolutions, the Ohmic heating changes the interior 
structure and thus decreases r\. As a result of the feedback, the skin depth becomes smaller and thus the 
induced electromagnetic effect is mostly confined in the radiative layer. 
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Table [I] shows that the stellar obliquity A remains zero in Cases 1-3 as expected, because 
(T z i) is the only component of the dissipation torque acting on the stellar spin. Moreover, A 
in Cases 4 and 6 change similarly; AA ~ —11° for both cases, meaning that the dissipation 
torques turn the system toward the spin-orbit alignment in Case 4 and toward the anti- 
alignment in Case 6 at similar rates. By contrast, A hardly alters in Case 5 when A» = 90°. 
In the presence of a, the evolution of A is more complicated; the dissipation torques can either 
excite or damp the stellar obliquity. Table [1] shows that A in Cases 7 and 8 is decreased by 
about 2°-3°, while A in Case 9 is increased by about 3°. The changes of A are nonetheless 
much slower than those in Cases 4 and 6. 

The aforementioned evolutions of A for Cases 4-9 can be understood through exami- 
nation of Equation (IT7j) . Because a = in Cases 4-6, dL spin /dt = and dL orb /dt < in 
these casea^ as has been shown in §2. Equation (fT7j) then indicates that the dissipation 
torques acting on the star are unable to spin up/down the star but contribute entirely to the 
evolution of A. This explains why the evolutions of A in Cases 4 and 6 are faster than those 
in Cases 7, 8, and 9. Besides, we find that the orbital axis moves faster toward the stellar 
spin instead of the other way around. It is due to the fact that M p D 2 n/ I^u* ~ 0.1-0.01 <C 1 
in our model. Namely, it is the terms associated with l/(M p D 2 n) rather than with 1/(7*0;*) 
on the right-hand side of Equation ffTTl) (or equivalently Equation ffl5l) ) that dominate A. As 
a result, Equation ffTTl) gives similar decreases in A in Cases 4 and 6. On the other hand in 
Case 5, Aj = 90° and hence only the terms with l/(i*u;*) in Equation ( fl7|) exist at the be- 
ginning, which is negative and small. This explains why A hardly evolves in Case 5; namely, 
AA is only —0.02°. We can set Equation f[TT|) equal to zero in the case of a = and find the 
equilibrium orientation for the stellar spin, which gives A = arccos(— M p D 2 n/I*u*) ~ ±90°, 
consistent with the small A in Case 5. The equilibrium is unstable as we can anticipate 
from Cases 4 and 6. Roughly speaking, the dissipation torque turns the stellar spin and the 
orbital axis toward alignment when |A| < 90° and toward anti-alignment when > 90°. 

When a ^ 0, dL spin /dt becomes non-zero. In the case of n > a;*, dL spin /dt > and 
dL or b/dt < 0. Equation ([17]) implies that for Cases 7, 8, and 9, the term associated with 
dL orb /dt damps A, while the term with dL spin /dt excites A. When a = 0, dL spin /dt = 
and therefore the excitation term vanishes. As a increases from zero, the excitation term 
increases as well. This reiterates the point described in the above paragraph that A damps 



5 dL sp i n /dt — here should be distinguished from the assumption cj* = constant. dL sp i n /dt = in 
Equation (Tl7|) results from internal torques in the diamagnetic interaction between the star and planet for 
the cases of a = 0. On the other hand, the assumption that w* = constant is made under the consideration 
of any external angular momentum transfer between the T Tauri star and the environment, such as via disk 
locking. 
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much more slowly when a = 45°. When a = 90°, it turns out that the excitation term 
dominates over the damping term, leading to A > in Case 9. 

As we have seen, the values of |A| in Cases 4 and 6 change similarly. Indeed, inspection 
of Figure [7] reveals that the evolutions of R p and D in Cases 4 and 6 are very similar. As 
has been described above, the reason is that their secular interactions are almost the same 
when the stellar dipole is aligned or anti-aligned with the stellar spin (i.e. a = 0). 

Cases 10-14 present the evolutions of the planet placed initially at the farther distance 
Di pa 0.03 AU from its T Tauri star. Hence, the magnetic interaction and the resulting 
Ohmic dissipation are weaker at the beginning of the evolution (i.e. Qohmic ~ 10 29 erg/s) 
than those in Cases 1-9. Consequently, the planet's orbit decays to D < 0.018 AU over a 
timescale of a few Myrs as shown in Figure El Unlike the initial rise of R p in Cases 1-9, 
R p first decreases quickly over a short period of time, ~ 10 5 years in these cases. Then the 
planet undergoes a slow change of the size over a few Myrs, meaning that Qohmic ~ L in 
Equation ([6]) such that the radius is more or less maintained by the Ohmic dissipation at 
R p ~ 2-2.04i?/. In addition, Figure |S] and Tabled] show that D, R p , and A of Cases 11 & 12 
evolve almost identically, as expected for the same reason as the similarity between Cases 
4 & 6. The stellar obliquities in Cases 11 & 12 decrease by a similar amount, about 20° at 
the end of the simulation, but decrease only by about 5° in Cases 13 & 14. The results are 
consistent with the comparative study for Cases 4, 6, 7, and 8; i.e., A decays faster when 
a = in these cases. The faster decreases of A in Cases 11 & 12 give rise to the smaller 
heating rates than those in other cases during the late stages of the evolution, leading to 
faster contraction and slower orbital decay after t pa 6-7 Myrs. As a result, the planet in 
Cases 11 and 12 is unable to reach Roche-lobe overflow at the end of the simulation, whereas 
the planet in Cases 10, 13, and 14 migrates in fast enough to finally fill the Roche lobe during 
the T Tauri phase. 

In Cases 15-19, we consider a young hot Jupiter initially on the even larger orbit at Di pa 
0.04 AU. These cases correspond to the scenario that the planet lies in a large magnetospheric 
cavity. The larger the orbital distance is, the weaker the stellar magnetic fields and the 
interactions are. In all cases, the Ohmic dissipation rates are around 10 28 erg/s throughout 
the simulation. As a result, the young planet first contracts and then gradually attains 
quasi-thermal equilibrium at small radii R p pa l.QRj after t ~ 4-5 Myrs, as illustrated in 
Figure |9j The orbits do not decay to the distance D < 0.039 AU from the star. Hence, 
the planet never reaches its Roche radius at the end of the simulation. The corresponding 
changes in A shown in Table [TJ are much smaller as well in these cases. 



Laine et all (120081 ) focused on a planet that already fills its Roche lobe at Di = 0.04 AU 
(i.e. R p pa 5Rj), and calculated the resulting dissipation rate and mass loss rate. The interior 
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structure of the planet is composed of an isothermal envelope and a polytropic core in their 
study, which does not truly take into account the energy equation. Although they suggested 
that the dissipation can inflate the planet and trigger mass loss through Roche lobe overflow, 
the authors cautioned that the initial condition of the Roche-lobe filling plane t would not 



be val id. In Case 20, we apply our interior structure to the fiducial model of lLaine et al. 
( 120081 ) with the initial R p = 2.045i?j at t = ti. The evolutions of R p and D are plotted in 
Figure [10J Since the planet is located outside the co-rotation orbital radius, the dissipation 
torque induced by the forcing w_ > drives the planet to migrate outwards. Although the 
planet in this case is less massive (M p = 0.63Mj) than in other cases (M p = lMj), the 
Ohmic heating i s not strong enough to inflate the less massive planet to its Roche radius 
as suggested by lLaine et all ( 120081 ). Rather, the planet radius remains almost constant 
throughout the T Tauri phase. 



5. Summary and discussions 

We revisit the magnet ic interaction between a hot Jupiter and its T Tauri star investi- 



gated by lLaine et al.l ( 120081 ). In the original work, the authors considered a Roche- lobe sized 
hot Jupiter without its own magnetic fields. The stellar spin was assumed to be aligned 
with the orbital axis; i.e. the stellar obliquity A = 0. As the planet orbits its parent star, 
the Ohmic dissipation in the planet is induced by the stellar magnetic dipole tilted away 
from the stellar spin with an angle a. To calculate the electric resistivity of the planet, the 
interior structure was modelled as a sphere consisting of a polytropic core and an isothermal 
outer layer. In their model, the planet lies outside the co-rotating orbital radius such that 
the forcing frequency w_ > 0. Based on their fiducial model, the authors suggested that 
the dissipation torques are not able to cause any significant orbital change. Nevertheless, 
the Ohmic heating occurring in the outer part of the planet is intense enough to inflate the 
planet up to the Roche radius. The mass loss through the Lagrangian 1 point toward the 
central star provides the angular momentum to the planet and thus possibly halts the planet 
migration in the disk. 

Motivated by a wide range of the stellar obliquity detected in hot- Jupiter sys tems (e.g . 



Winn et al. 



2011). which in theory could be excited during the T Tauri phase ( jLail 12012 



Foucart k, Laill201ll ). we extend the original model by considering the coupled evolution 
of the interior structure, planet's orbit, and the stellar spin in the presence of the stellar 
obliquity A. We focus on the secular evolution due to the dissipation torques and show that 
the Ohmic dissipation in the planet can be contributed linearly from the forcing associated 
with 4 frequencies: u + , a;*, n, and Owing to the complication of the problem involving 
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multiple frequencies, we begin with a couple of test runs based on a given interior structure 
for the dissipation calculation, which are further validated by the general torque-dissipation 
relation given by Equation f[2"0"j) as well as the skin-depth estimation using Equation (jSJ). 

The coupled evolutions are then carried out for a number of cases listed in Table [1] for 
the purpose of parameter studies. The evolutions are computed from t = 0.7 to 10 Myrs. 
The radius of the coreless hot Jupiter of lMj is 2.045i?j at the beginning. Initially, the T 
Tauri star is assumed to spin at a rate uo* = n/2 to imitate the final stage of the planet 
migration scenario with a giant planet inside the magnetospheric cavity of the disk (Lin et 
al. 1996; Rice et al. 2008; cf. Bemtez-Llambay et al 2011). Since the planet lies inside the 
co-rotating orbit, the planet continues to migrate inwards due to the Ohmic dissipation in 
the planet. Without modelling the star-disk magnetic interactions, the co-rotating orbit is 
simply assumed to correspond to the inner edge of the magnetically truncated disk despite 
the presence of A and a. This is certainly one of the limitations of the study. In most of the 
cases, uj* is assumed constant to simply resemble any processes, such as disk locking, that 
maintain the stellar spin. 

Three initial orbital distances Di rj 0.02, 0.03, and 0.04 AU are considered. With our 
input parameters for Di w 0.02 AU, the intense dissipation confined near the planet surface 
only enlarges R p by < 3%. Nonetheless, the dissipation torques decay the orbit to the Roche 
zone in just a few 10 5 to about 1 million years. The torques also evolve A when \ ^ 0. Since 
M p D 2 n <C A is primarily contributed from the movement of the orbital axis rather than 
the stellar spin axis. When a is zero, there exists an unstable equilibrium for the orientation 
of the stellar spin axis, which points roughly about 90° from the orbital axis. Consequently, 
the dissipation torques direct the orbital axis toward the stellar spin for a prograde orbit but 
toward the anti-parallel direction to the spin for a retrograde orbit. When a is non-zero, the 
orbital axis and the stellar spin can either evolve toward alignment/anti-alignment for small 
a or become more misaligned for a < 90°. Because the stellar spin is not spun up/down by 
the dissipation torque when a = 0, it follows that the stellar obliquity evolves more quickly 
when the stellar spin is parallel/anti-parallel to the stellar magnetic dipole. 

For the young hot Jupiter initially at the farther distance Di w 0.03 AU, the dissipation 
is modest but still strong enough to more or less sustain the initial R p except for the cases 
with a = 0. The relatively fast decrease of A for a = weakens the dissipation and the 
resulting torques, leading to the planet contraction and slow orbital decay in the late stage 
of the evolution. Therefore, in terms of the cases we have studied, the planet with a ^ 
undergoes substantial orbital decay in a few Myrs and finally overflows the Roche-lobe, 
while the planet with a = can shrink its orbit but not sufficiently to allow for Roche-lobe 
overflow. 
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Owing to the weaker interaction at the larger orbital distance, the planet of 1 Mj in 
all the cases starting from Di w 0.04 AU contracts and then roughly reaches quasi-thermal 
equilibrium during the T Tauri phase, with the final R p smaller than those in the cases for 
Di 0.02 and 0.03 AU. The corresponding orbital decays and the changes in the stellar 
obliquity are substantially smaller. The planet moves barely from its initial orbit and thus 
is not able to r each its Roche-lobe. We also carry out the simulation for the fiducial model 



in 



Laine et al.l (120081 ) and find that the Ohmic heating can only sustain the radius of the less 
massive young hot Jupiter (M p = 0.63Mj), rather than thermally expanding the planet to 
its Roche radius as suggested by the authors. 

The induced dissipation rates are as high as 10 30 ~ 31 erg/s when the planet moves to 
about D < 0.02 AU. The intense heating near the planet's surface does not generate local 
temperat u re inv ersio n in our model, i n con trast to the surface heating models presented in 
Gu et al.l (120041 ) and lWu fc Lithwickl (120121 ). It is probably because the dissipation respon- 
sible for the temperature inversion lies within a thin shell fairly deep in Gu et al. 2004 
(a prescribed narrow Gaussian region) and in Wu & Lithwick 2012 (a region at about the 
optical depth of 100), while the dissipation in our diamagnetic induction model is deposited 
so close to the surface that it is easily lost and so there is no local maximum in T. 

In this work, we introduce the planet at £j = 0.7 Myrs and adopt a constant magnetic 
dipole moment. In theory, a gas giant planet can form later and migrate to the mag- 



netospheric cavity at a later time (e.g., Ilda &: Linl 120081 ; iMordasini et a. 



200 91). Besides, 



the magnetic dipole mom ent may decay over the course of a few Myrs (jJohns-KrullI 120071 ; 
Yang fc Johns-Krulll 120111 ). Hence, our results probably give the suggestive values for the 
maximum changes of the stellar obliquity, orbital distance, and planet radius during the T 
Tauri star phase. 

The orbital decay of a young hot Jupiter in our magnetic model is not significant unless 
the distance to the T Tauri star is sma ller than about 0.03 AU. Some other process, such as 
gravitatio nal tides (IChang et al.l 120101) or a small magnetospheric cavity during FU-Orionis 



outbursts (IBaraffe et al 



2009 



Adams et al.ll2009i ). can bring the planet in to that distance. 



The Ohmic mechanism could be the final stage in bringing a planet in very close to the T 
Tauri star or even leading to a Roche-lobe overflow. This may provide one of the explanations 
for the pile-up of hot Jupiters with the orbital periods of ~ 3 days, and could also reduce the 
too-h igh population of hot Jupiters inferred from the population synthesis model (llda fc Lin 
20081 ). In this work, we do not study the post-evolution of a Roche-lobe filled planet. In 
terms of our model parameters, a young hot Jupiter located within 0.03 AU from its T Tauri 
star can undergo fast orbital decay on the timescales much shorter than 10 million years. 
Therefore it is possible that the young planet overflows its Roche lobe, migrates out, then 
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migrates in, and overflows again. Consequently, the planet suffers from intermittent mass 
losse s until its density i s low enough to go through the stage of the runaway adiabatic mass 
loss (jChang et al.ll2010l ). leading to the demise of the planet during the T Tauri phase. 



As has been described in the Introduction, there is a large body of literature devoted to 
a variety of magnetic interactions between a hot Jupiter and its parent star, some of which 
can also cause the angular momentum to transfer between the planet's orbit and the stellar 
spin. The efficiency of angular momentum transfer is model dependent, relying on the Ohmic 
dissipation rate. In the T Tauri phase, the presence of a disk i s exp ected to magnetically 
affect the stellar spin and perhaps the planet's orbit. lLail ( 120121 ) and iFoucart &: Lail (1201 if ) 
considered a hybrid magnetic model including diamagnetic induction and magnetic-field 
linkage for the purpose of the generati o n of the stellar obliquity. In the study presented 
here, we follow the work by lLaine et al.l (120081 ) and therefore focus only on the diamagnetic 
interaction between the planet and its T Tauri star. Our simple model suggests that the 
stellar obliquity starting from a non-zero value may further evolve after the planet migrates 
into the magnetospheric cavity, making the orbit of the young hot Jupiter incline with the 
disk plane. As a result, a hot Jupiter does not necessarily lie on the same orbital plane with 
the planets farther out from the central star. Whether or not our model can provide a wide 
range of stellar obliquities at the end of T Tauri phase depends on the initial distribution of 
stellar obliquity as well as the distribution of the direction of stellar dipole moment relative 
to the stellar spin. 

In the Introduction, we also caution that the skin depth beneath the photosphere is 
one of the major uncertainties of the model. In the presence of a planetary ionosphere or 
magnetosphere, the value may be appreciably smaller than what w e comput e in this work. 
Nevertheless, as an analog of the star-disk magnetic interactions (jLail 119991 ) . our diamag- 
netic model and other magnetic interactions should be considered together for the orbital 
evolution inside the magnetospheric cavity. While theoretical models are under develop- 
men t, it is conceivable in the future that photometric variability on timescales of a few days 



2008 



(e.g. Bouvier et al.ll2007l) and spectropolarimetry applied to T Tauri stars (e.g. iDonati et al. 



Long et al.ll201lf ) would serve as possible detection methods, to search for the variabil- 



ity modes and magnetic perturbations that are associated with the orbital motion of such a 
young hot Jupiter during the T Tauri star stage. 
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A. magnetic scalar potential in the presence of stellar obliquity 

We work on the problem with two sets of the coordinate sy stems O'x'y'z' and Oxyz 



illustrated in Figure (TJ which were employed in lLaine et al.l (120081 ) with the orbital angular 
velocity and stellar spin around the z (and z' axis). The magnetic dipole moment m is tilted 
with the angle a from the stellar spin and thus can be expressed as m = (/z /47r)(I?*ii^/2)rh 
with the unit dipole moment vector rh = sin a cos u)*ti' + sin a sin + cos ak! as viewed in 
the star's frame. Here B* is the surface stellar field and R* is the stellar radius. Then let the 
stellar spin rotate along the x' axis clockwise (counter-clockwise) on the y'-z' plane such that 
the stellar obliquity angle is A > (< 0). This gives u>* in the inertial frame of our problem. 
Note that the sign definition of A agrees with that used for the Rossiter-McLaughlin effect. 

We then rotate the stellar spin along the z' axis clockwise at the angular velocity n to 
give the stellar spin as observed in the co- moving frame of the planet. As a result, the unit 
vector of the stellar magnetic dipole moment rh ob i iquity as viewed in the planet's rest frame 
becomes 

(sin a (cos uu*t cos nt + cos A sin uij, sin nt) + sin A cos a sin nt 
sina(— cosco>*£ sinnt + cos A sin u;*i cosnt) + sin A cos a cos 
— sin A sin a sin to J, + cos A cos a 

(Al) 

The rotation matrices involved in the above equations are 

\ / cosnt sinnt 

Pv' = I cos A sin A , P z ' =1 — sin nt cos nt 
- sin A cos A / \ 1 



Note that when the obliquity A = in Equation flAip . we recover rh = sin a cos wti' + 
sin a sin out}' + cos ak with the Doppler-shifted frequency u> being u, - n as viewed in the 
co- moving frame of the planet. Moreover, if the obliquity is retained but the magnetic axis is 
aligned with the stellar spin (i.e. a = 0), we obtain m ob ii quity = sin A sinnti' + sin Acosntj' + 
cos Ak as should be expected in the co- moving frame of the planet. 

Then we can express the magnetic scalar potential due to the stellar magnetic dipole 
moment with non-zero stellar obliquity in the co-moving frame of the plane as follows (cf. 
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f- p _ m |-f2 Le sm2<^sina((cosa; + t + cosu;_t) + cos A(cosw_t — coscu+t)) 



Campbell! Il983[ ): 

V ^° m v' rh 

4-Kr' 3 obhquity 

= r sin 9 sm <p> sin afcos u;*r cos nt + cos A sin nt sm u;*£) + r sin sin 09 sin nt sin A cos a 

Anr IA 

+ (D — r sin # cos 99) sina(— sin rat cosu^t + cos A cos nt sinu^t) 

+ (D — r sin 9 cos (p) sin A cos a cos nt + r cos #(— sin A sin a sin u^t + cos A cos a)] 

HqUIT i 1.. 

~ - — "i - sm^sma((cosw + t + cosa;_t) + cos A(cosu;_t — cosw+t)) 
+ sin ip sin rat sin A cos a 

+ cosy? sina((sinu;_t — sinw+t) + cos A(sina; + t + sinu;_t)) 

+2 cos sin A cos a cos nt] — P® sin A sin a sin u;*i} 

3u mr 2 2 1 

8ttP 4 1 2[ 6 

H — sin 2 cp sin nt sin A cos a 
3 y 

+ - cos2o9sina((sinw_t — sinw + t) + cos A(sinw + t + sinw_t)) 
+— cos lip sin A cos a cos nt] 

— P® [— sina((sinw_t — sinw+t) + cos A (sin u + t + sinw_t)) + sin A cos a cos rat] 
2 

+-P2 1 cos sin A sin a sinc^t} (A2) 

where cj + = + ra, cj_ = — n, and the approximation of V* is obtained by the expansion 
of r' -3 up to the order of (r/D) 2 . In the above equation, the time-independent terms are 
dropped out because they do not contribute to the electromagnetic induction. Besides, Pj 
are associated Legendre functions, defined bj@ P° = cos 9, Pi = — sin 9, P 2 ° = (1/2) (3 cos 2 9— 
1), p 2 i = -3 sin^ cos 9, and P| = 3sin 2 fl. 



6 Note that one can obtain the expression for m < from m > using the relation P ; m (a;) = 
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Using V* = —d r <f)*, we have the poloidal scalar of the stellar field as follows 



2 ,irl 



^° „ {-Pito sin sin a((cosu + t + cos w_t) + cos A(cosw_t — cosw + t)) 
+ sin (p sin rat sin A cos a 

+ cos sin a((sinu;_t — sinw+t) + cos A(sinu; + t + sincj_t)) 
+2 cos sin A cos a cos rat] + P° sin A sin a sin u;*t} 



H^mr o 1 

+ — — " 2 [— sin 2<£> sin a((cosa; + t + cosu;_t) + cos A(cosu;_t — cos u; + t)) 



6 



+— sin 2ip sin rat sin A cos a 

+ - cos 2ip sin a((sinw_t — sin ui+t) + cos A(sinw + t + sinw_t)) 
1 

+- cos 2<£> sin A cos a cos rat] 

+P® [- sin a((sinu;_t — sinw+t) + cosA(sinu; + t + sinu;_t)) + sin A cos a; cos rat] 



2 

-P\ cosy? sin A sin a sin u^t}. 
3 



(A3) 



Note that we recover the original form of the poloidal scalar 0* in lLaine et al.l (120081 ) when 
A = 0. In addition, because there is a poloidal scalar outside the planet (fi p generated by the 
fields induced by 0* inside the planet, P has the same time and angular dependence as (ft*. 
Hence, P is given by 



(Si sin u*t + <5 2 cos u*t) 



<p p (r > Rp,t) 

1 

^-wsino? . . . 

+/xo"i [ i a i sm u+t + a 2 cos u + t + a 3 sm w_ + a 4 cos w_t + a 5 sin rat + a 6 cos nt) 



COS if 



r 



a-j sin u + t + a 8 cos u + t + a 9 sin U-t + qji cos U-t + an sin nt + ai2 cos rat)] 



U P° 

H sin w + t + /3 2 cos tu + t + /3 3 sin uj_t + /3 4 cos w_t + /3 5 sin rat + /3 6 cos rat) 

1 COS if 

+/i -P 2 2~ ^ S111 + £ 2 COS 

2r sin2<z>. 

+/io° 2 [ — 2 — w 1 sm UJ +t + 7 2 cos u +t + 73 sm + 74 cos u -t + 75 sm n t + 76 cos ^) 



cos 2ip 



(77 sin w+t + 7s cos u + t + 79 sin uj_t + 7 10 cos w_t + 711 sin rat + 7 12 cos rat)](A4) 



The above form for 0* along with the separation of variables for the solutions to the 
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induction equation suggest that the induced potential scalar <p p can be in general given by 



iuj—t 



(A5) 



When A = 0, the above expression is reduced to only one term that is associated with w_. 
Inside the planet (r < Rp), the r-dependence of (fi p is given by 

r d 2 



<ir 2 



Z(Z + 1) zu; 
r 2 



G {)W (r) = 0, 



(A6) 



where the forcing frequency u denotes co + , w_, w*, or n. If 77(7*) is known, the above equation 
with each frequency can b e solved by rearrang ing the equation to 4 first-order ODEs as have 
shown in Equation (10) of lLaine et al.l (120081 ). subject to the boundary conditions 



G't(Rp 



I 



Gi(Rp)-(2l + l)R l p = 0, 



G[{r « 0) - l -^Gi{r w 0) = 0, 



(A7) 
(A8) 



After obtaining the solutions for Gi }U) , we are ready to solve for the coefficients in the 
expressions of 4> p {r < R p ,t) and <j) p (r > R p ,t) by demanding the condition that the potential 
scalars and their derivatives should continue at r = R p for each frequency. Namely, + 



■ > p,out 



J p,in 



and d r 



:+< 



'p,out ) 



0, 



rfp,in 



air = R p for each frequency. As a result, the terms 



associated with u + contribute 10 algebraic equations from Y^ l e luJ+t , Y^e^*, and Y^ 2 e tU)+t . 
Likewise, the terms depending on w_ and n give rise to 10 algebraic equations each. On 
the other hand, the terms with only contribute 6 algebraic equations from Y^e w * t and 



y 2 ±1 e ia; * < , but keep in mind that the two terms associated with siny? from Y^ 1 ^* 1 should 
vanish as they do not exist in <p p , leading to the relation between C\ u and C^ 1 w (i.e. 

w = C^ 1 ) and hence reducing to 4 equations. Including their derivatives counterparts, 
we shall solve 20 algebraic equations for the 20 coefficients associated with oj + , lo_, and 

Besides, 8 equations for the 8 coefficients associated with cu*. In the end, we have 68 
algebraic equations at r = R p and solve for the 68 coefficients, which are 5\ t 2, ai,---i2, Pi,-6j 



e li2 , 71,-12, Re(Cf ^_ >n ), Im(Cf ^_ >n ), Re(C° ^_ >n ), bn(^^_,„). *MCf „ 
lm{Ci 



±i 



u>+,co— ,n/i 



Re(C? )W J, Im(C? J, Re(C*J, and Im(C^J. Then we know p , out and 



J p,m- 



7 T he 20 equations for ui_ are equivalent to the set of 20 linear equations in the Appendix B of lLaine et al. 
(|2008l ). There are typos on the right-hand sides of their 5th to 7th linear equations. 
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B. Secular evolutions of the precession and obliquity of the stellar spin 



While the stellar spin and the planet's orbit exchanges angular momentum due to the 
dissipative torques in our magnetic model, the total angular momentum vector is conserved. 
Bearing this in mind, we have the secular evolution of the precession angle governed by (cf. 
Goldreich fc Pealelll97oh 



sin X^I^u* 



dt 



(T x "), 



(bi; 



where A* is the angle between the stellar spin and the total angular momentum. Since the 
stellar spin angular momentum, orbital angular momentum, and total angular momentum 
form a vector triangle, it is straightforward to show from the triangle that 

sin A 



sin A, 



^/{hujM p D 2 nf + 1 + 2(hoojM p D 2 n) cos A 
This equation yields the expression of /(A) in Equation (JT31 . 



= /(A). 



(B2) 



Now we turn to the evolution of the stellar obliquity. We define d\*/dt and d\ n /dt as 
the contributions to the secular evolution of A due respectively to the stellar spin and orbital 
angular momentum moving toward/away from the total angular momentum. As has been 
described in the main text, A is caused by the components of (T y //) and (T z i) normal to the 
stellar spin and by the back reaction —(T v >>) ac ting to the orbital angular momentum. Hence 
we have (cf. iGoldreich fc PealeHl97(j Ea7ll999h 

(Tyii) cos A — (T z i) sin A, (B3) 



dt 



Therefore the equation 



71 f n 2 d\ n 
M V D n—— 

p dt 

G?A* 



dX 
~dl 



dt 



dXg 

dt ' 



(B4) 



(B5) 



gives the expression in Equation (1151) . 



When M p D 2 n > /.w*, A, 



A and A* ~ A as they ought to be because the total 
angular momentum is almost contributed from the orbital angular momentum. However, 
M p D 2 n < I*u* for a system consisting of a hot Jupiter and a T Tauri star. As a result, the 
planet's orbit can evolve more significantly than the stellar spin. 
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Table 1: model calculations, m = 4 x 10 34 A m 2 , M p = lMj, M* = 1M , L* = L Q . The 
subscript i denotes the initial value, while the subscript / means the final value due to either 
the planet reaching its Roche lobe or t = t end . In addition, at t iy R p = 2.045-Rj, (0") = 0. 
The column "overflow" indicates in which cases the planet reaches the Roche-lobe overflow 
before t = t en d- The figure number for the results of each case is indicated in the column 
"Figure" 



Case 


A (AU) 


a 


Aj 


A/ 


overflow 


Figure 


1 


0.02 


45° 


0° 


0° 


yes 


m 


2 


0.02 


45° 


0° 


0° 


yes 


7 


3 


0.02 


90° 


0° 


0° 




7 


4 


0.02 


0° 


45° 


33.77° 


yes 


7 


5 


0.02 


0° 


90° 


89.98° 


yes 


7 


6 


0.02 


0° 


-135° 


-146.02° 


yes 


7 


7 


0.02 


45° 


45° 


43.26° 


yes 


7 


8 


0.02 


45° 


-135° 


-137.97° 


yes 


7 


9 


0.02 


90° 


45° 


48.25° 


yes 


7 


10 


0.03 


45° 


0° 


0° 


yes 


8 


11 


0.03 


0° 


45° 


21.97° 


no 


8 


12 


0.03 


0° 


-135° 


-157.89° 


no 


8 


13 


0.03 


45° 


45° 


40.56° 


yes 


8 


14 


0.03 


45° 


-135° 


-141.13° 


yes 


8 


15 


0.04 


45° 


0° 


0° 


no 


9 


16 


0.04 


0° 


45° 


44.44° 


no 


9 


17 


0.04 


0° 


-135° 


-135.55° 


no 


9 


18 


0.04 


45° 


45° 


44.88° 


no 


9 


19 


0.04 


45° 


-135° 


-135.16° 


no 


i 


20 a 


0.04 


90° 


0° 


0° 


no 


ma 



"Fiducial model in lLaine et al 



See the text for the details. 
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Fig. 1. — The coordinate systems adopted in this study for calculations. The Oxyz coordi- 
nate system is anchored at the planet, with the origin O at the planet's center and x-axis 
pointing to the proto-star. On the other hand, the origin O' of the O'x'y'z' system lies at the 
proto-star's center with the y'-axis pointing to O. The coordinate systems are the same as 
those in Laine et al. (2008). At t = the stellar spin uj* is placed on the y'-z' plane and is 
inclined at the obliquity angle A relative to the direction of the orbital axis n (i.e. the vertical 
axis z'). As is viewed in the Oxyz system co-moving with the planet, u;* rotates around the 
vertical axis with the angle A at the rate of the orbital angular frequency n. Moreover, as 
the star spins, the stellar dipole moment m rotates around u;* with the misaligned angle a. 
We let the y' axis always secularly follow the precession of uj* such that u>* always lies on 
the y'-z' during the secular evolution. 
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Fig. 2. — Dependence of the Ohmic heating rate on A, a, and u;*/n. The Ohmic heating 
rates are calculated from Equation([7]). The mass, radius, and orbital radius of the young 
hot Jupiter are 1 Jupiter mass, 1.84 Jupiter radii, and 0.02 AU for this set of calculations. 
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Fig. 3. — Ohmic heating rate profiles in a young hot Jupiter of 1 Jupiter mass and 1.84 
Jupiter radii, a = 10° and 90° for the stellar obliquity A = 45° are plotted. As illustrated 
in the figure, most of the Ohmic dissipation occurs in the outer part of the planet because 
the induced electromagnetic fields can only penetrate from the surface down over a length 
scale comparable to the skin depth ~ ■^JrjJuJ. The radiative-convective interface is located at 
m 1.295 x 10 10 cm. Note that the total heating profile overlaps with the profiles associated 
with other forcing frequencies (see the text in §3 for the details). 
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Fig. 4. — Same as Figure [3] but for a planet in a retrograde orbit with A = 135°. 
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Fig. 5. — Same as Figure [2] but the Ohmic heating rates are computed from Equation ( l20l) . 
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Fig. 6. — Comparison between the total heating profile derived from Equation (J2J) and that 
from the full version of the Saha equation for the cases corresponding to the top left (left in 
this figure) and bottom right (right in the figure) panels of Figure [3j 
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Fig. 7. — Coupled evolutions of R p and D of the young hot Jupiter in Cases 1-9. The number 
labelled next to each curve is the case number. The evolutionary curves for Cases 4 and 6 
are very similar. 
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Fig. 8. — Coupled evolutions of R p and D of the young hot Jupiter in Cases 10-14. The 
numbers labelled next to the curves indicate the different cases. The evolutionary curves for 
Cases 11 and 12 almost overlap. 
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Coupled evolutions of R p and D of the young hot Jupiter in Cases 15-19. 
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Fig. 10. — Coupled evolutions of R p and D of the young hot Jupiter in Case 20. 



